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ABSTRACT  OF  THE  DISSERTATION 

Observation  of  High  Frequency  Internal  Waves 
in  the  Pacific  Ocean 

by 

Robert  Libman  Zalkan 

Doctor  of  Philosophy  in  Oceanography 
University  of  California,  San  Diego,  1968 
Professor  Fred  N.  Spless,  Chairman 

Observations  of  high  frequency  Internal  waves  were 
made  from  1  to  14  December  1966,  in  the  deep  sea  off  3a Ja 
California.   The  fluctuations  of  the  depth  of  an  isotherm 
were  measured  with  a  three-element  horizontal  array 
attached  to  a  stable  platform. 

The  waves  are  characterized  as  a  broad  band 

viii 


phenomenon  with  a  continuous  distribution  in  frequency.   For 
short  Intervals  of  time,  a  narrow  frequency  band  within  the 
continuum  is  adequately  described  as  a  horizontally  plane 
wave  of  a  single  vertical  mode.   Furthermore,  this  simpli- 
fied structure  is  stationary  over  time  spans  of  several 
days.   High  modes  are  present  in  the  low  frequency  waves. 
Above  4  cycles  per  hour,  however,  the  first  mode  predomi- 
nates. 

The  spectral  shape  is  consistent  with  the  shear 
limited  equilibrium  spectrum  proposed  by  Phillips  (1966). 
In  addition,  the  predominance  of  the  first  mode  at  high 
frequencies  further  emphasizes  the  importance  of  shear 
instability  in  internal  wave  propagation. 

The  horizontal  properties  of  the  wave  field  indi- 
cate well-defined  directions  of  narrow-beamed  propagation. 
These  directions  and  the  dispersive  properties  of  the 
propagation  have  led  to  the  identification  of  local  topo- 
graphic features  as  generating  areas  of  internal  waves. 
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CHAPTER  I 
INTRODUCTION 

Since  internal  waves  constantly  change  the  interior 
structure  and  velocity  fields  of  the  ocean,  they  may  have 
Important  effects  on  dynamic  height  calculations,  current 
measurement,  mixing,  marine  organisms,  underwater  sound  pro- 
pagation, and  submarine  operations.   An  adequate  description 
of  the  wave  field  is  required  in  order  to  determine  the  ex- 
tent of  these  effects.   In  addition,  an  understanding  of  the 
mechanisms  of  generation  and  decay  is  necessary  to  the 
determination  of  the  role  internal  waves  play  in  energy 
transfer  processes.   Specifically,  the  description  should 
Include  modal  structure,  energy  spectra,  and  directional 
information. 

At  the  present  level  of  theoretical  development,  a 
model  that  adequately  describes  these  features  is  clearly 
not  feasible,  and  little  is  known  about  the  mechanisms  of 
generation  and  decay.   Our  purpose,  therefore,  is  to  seek 
such  a  description  from  observation  of  high  frequency 
internal  waves  in  the  deep  sea. 

The  complex  structure  of  internal  waves  makes  both 
data  collection  and  analysis  difficult.   Unlike  surface 
waves,  internal  waves  do  not  have  a  unique  dispersion 
relationship.   Instead,  there  is  a  different  relationship 


for  each  member  of  the  Infinite  set  of  possible  modes.   In 
addition,  the  maximum  vertical  displacements  of  these  waves 
occur  at  intermediate  depths  rather  than  at  the  free  sur- 
face, and  are,  therefore,  difficult  to  observe. 

To  overcome  these  difficulties,  simultaneous  measure- 
ments of  the  fluctuating  depth  of  an  isotherm  were  recorded 
with  a  three-element  horizontal  array.   This  procedure 
allows  the  establishment  of  a  phase  difference  (travel  time) 
relationship  leading  to  a  vector  wave  number  and,  therefore, 
directional  determination.   The  wave  number  magnitudes  can 
then  be  tested  against  the  theoretical  dispersion  curves  for 
confirmation  of  progressive  wave  phenomena. 

Stokes  (1847)  first  showed  that  Internal  waves  can 
exist  in  a  stratified  fluid  consisting  of  two  homogeneous 
layers.   Later  Webb  (1884)  and  Greenhill  (1887)  dealt  with 
the  multi-layered  case.   In  1933,  FJeldstad  extended  the 
theory  to  include  a  continuous  density  gradient.   More 
recent  systematic  treatments  of  the  subject  are  by  Eckart 
(I960)  and  Tolstoy  (1963).   A  summary  of  the  theory  of  small 
amplitude  Internal  waves  in  still  water  is  given  in 
Appendix  A. 

Several  Investigators  have  attempted  to  obtain  di- 
rectional information  by  towing  a  thermistor  chain  in  a 
sequence  of  traverses  of  different  direction.   This  pro- 
cedure allows  one  to  obtain  a  directional  description  of 
the  wave  field  by  interpretation  of  measured  Doppler  shifts. 


Using  this  method,  Metcalf  et  al  (1962)  reported  internal 
waves  of  17  minute  periods  traveling  south  from  the  equator. 
However,  LaFond  (1964)  in  the  Pacific  and  Charnock  (1965) 
in  the  Atlantic,  had  little  success. 

Ufford  (1947)  appears  to  be  the  first  to  make 
simultaneous  measurements  of  internal  waves  at  spatially 
separated  points.   His  first  method  was  to  make  repeated 
lowerings  of  bathythermographs  from  the  bow,  midsection, 
and  stern,  of  a  ship  located  over  the  San  Diego  trough,  in 
1200  meters  of  water.   The  results  were  inconclusive,  due 
to  the  small  phase  differences  and,  therefore,  large  error 
limits  in  the  determination  of  wave  number.  To  increase  the 
phase  differences,  Ufford  moved  to  shallow  water  (40  meters) 
and  employed  both  bathythermograph  lowerings  from  three 
anchored  ships  and  a  buoyed  array  of  three  thermocouples. 
Unfortunately,  each  observation  spanned  only  a  few  hours  of 
time,  and  he  used  only  those  portions  of  the  records 
(average  1  hour,  7  minutes)  in  which  a  single  frequency 
predominated. 

Further  work  in  shallow  water  (18  meters)  has 
been  accomplished  at  the  Navy  Electronics  Laboratory's 
Oceanographlc  Tower.   Here,  a  three-element  array  of  iso- 
therm followers  has  been  successfully  used  to  determine 
directional  information  (LaFond  1962a,  Cox  1962).   Neverthe- 
less, mid-oceanic  observations  are  still  required  since 
shallow  water  propagation  differs  considerably  from  that  in 


deep  water.   In  particular,  the  wave  form, and  therefore  the 
frequency  content,  is  affected  by  the  nearness  of  the  sea 
floor  (LaFond,  196lb)  and  the  propagation  direction  by  shore- 
ward refraction  (LaFond, 1962b).   Also,  there  is  probably  a 
preference  for  the  first  mode  at  all  frequencies,  since  the 
shallow  depth,  and  thus  narrow  waveguide,  means  large  shear 
for  modes  greater  than  one. 

We,  therefore,  seek  our  information  from  the  deep 
ocean.   The  difficulty  of  employing  the  array  technique  in 
deep  water  is  attested  to  by  the  previous  absence  of  such 
measurements. 


CHAPTER  II 
OBSERVATIONAL  EQUIPMENT 

A  three  element  horizontal  array  was  obtained  by- 
attaching  heavy  duty  booms  to  the  Research  Platform 
FLIP's  (Fisher  &  Spiess,  1963)  bow  or  upper  structure.   The 
attachment  was  done  by  hinges  enabling  the  booms  to  be 
properly  stowed  while  FLIP  was  in  the  horizontal  position. 
Cables  to  the  sea-sensing  units  were  fed  over  the  ends  of 
the  booms.   Upon  "flipping"  to  the  vertical  position,  the 
booms  were  rigged  out  to  form  an  equilateral  triangle  of 
30.5  meters  (100  feet)  per  side  (Figure  1).   Each  boom  had 
its  own  independent  sensing  system  consisting  of  1)  an 
Isotherm  follower,  2)  a  vertical  staff  of  four  variable 
frequency  thermometers,  and  3)  a  vibrotron  pressure  guage. 

The  isotherm  follower  is  controlled  by  a  thermistor 
in  a  sensing  unit.   A  servo  system  is  arranged  to  vary  the 
depth  of  the  sensing  unit  in  such  a  way  that  the  temperature 
measured  by  the  thermistor  remains  constant.   The  thermistor 
signal  is  developed  across  a  DC  bridge*   One  leg  of  the 
bridge  is  connected  by  a  cable  to  a  thermistor  in  the  sea- 
sensing  element.   The  other  leg  is  composed  of  a  variable 
bank  of  resistors  for  selection  of  the  desired  temperature. 
When  the  thermistor  has  the  same  resistance  as  the  variable 
resistors,  the  bridge  is  in  balance.   When  the  temperature 
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at  the  thermistor  changes,  the  thermistor  changes  resistance 
and  a  DC  voltage  occurs  across  the  bridge.   This  voltage  is 
fed  through  magnetic  amplifiers  and  silicon  controlled 
rectifiers  to  the  windings  of  a  DC  motor.   The  polarity  and 
amplitude  of  the  signal  determine  the  direction  and  speed  at 
which  the  motor  turns.   The  motor  drives  a  winch  which 
raises  or  lowers  the  sea-sensing  element  until  the  bridge 
is  again  in  balance.   In  this  way  the  thermistor  in  the  sea- 
sensing  unit  is  kept  at  the  depth  of  an  Isotherm  and  hence 
the  name  Isotherm  Follower  (LaFond,  1961a;  Good,  1968). 

The  upper  part  of  the  sea-sensing  element  is  a 
hollow  pressure  case.  A  vibrotron  pressure  guage  was  instal- 
led in  this  portion  and  exposed  to  sea  pressure  via  an  oil- 
filled  chamber  with  a  flexible  diaphragm.  The  output  of  the 
vibrotron  was  used  for  determining  the  depth  of  the  sea- 
sensing  unit. 

Since  servo-systems  are  subject  to  dead  bands  (i.e., 
a  certain  level  of  signal  is  necessary  before  the  system 
responds),  a  staff  with  four  variable  frequency  thermometers 
(Squier,  1967)  was  installed  on  each  sea-sensing  unit  to  be 
used  as  a  vernier.    The  thermometers  were  placed  .6  meters 
(2  feet)  and  1.8  meters  (6  feet)  above  and  below  the  servo- 
system  thermistor.   A  scheme  of  interpolation  was  then  used 
to  determine  the  actual  depth  of  the  isotherm  (see  Appendix 
B). 

The  data  were  recorded  both  digitally  for  later 
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analysis  and  by  analog  methods  for  on  the  spot  monitoring. 
The  digital  data  collecting  system  consisted  of  1)  a  clock, 
2)  a  sequential  sampling  unit,  3)  a  counting  system,  and  4) 
an  Incremental  tape  recorder.   Every  twenty  seconds  the 
clock  put  out  a  pulse  which  triggered  the  sampling  unit. 
The  sampling  unit  consisted  of  a  stepping  switch  which  se- 
quentially connected  each  of  the  channels  of  Information  to 
the  counting  system  for  one  second's  duration.   The  counting 
system  was  a  period  counter  set  to  determine  the  number  of 
microseconds  required  for  4-096  cycles  of  signal.   This 
number,  along  with  a  serial  time  count,  was  then  recorded  on 
magnetic  tape  by  an  incremental  recorder.   The  information 
recorded  from  each  of  the  sea-sensing  units  consisted  of 
the  vibrotron  pressure  signal,  the  four  variable  frequency 
thermometer  temperatures,  and  the  vibrotron  pressure  signal 
again.   Each  follower  was  sampled  in  turn  and  then  the 
azlmuthal  orientation  of  the  array  was  recorded  to  conclude 
the  cycle.   After  the  19  pieces  of  information  had  been 
recorded,  there  was  one  second  of  dead  time.   Then  the 
clock  sent  out  a  pulse  and  the  whole  procedure  started  again. 

DC  voltages  for  analog  recording  were  obtained  by 
sending  the  three  vibrotron  signals  and  one  temperature 
signal  to  frequency  meters.   The  output  voltage  from  the 
meters  drove  the  pen  displacement  system  of  a  four  channel 
chart  recorder. 

Since  FLIP  is  free  to  rotate  about  the  vertical 


axis,  it  was  necessary  to  employ  an  azimuthal  orientation 
system.   An  error  signal  between  the  desired  heading  and 
FLIP's  actual  heading  and  a  rate  of  rotation  signal  from  the 
gyrocompass  repeater  were  used  to  drive  the  tilt  boxes  of 
two  hydraulic  motors.   The  motors  in  turn  drove  two  off-axis 
propellers.   The  torque  from  the  propellers  would  rotate 
FLIP  until  the  desired  heading  was  achieved.   If  FLIP 
drifted  off  the  desired  heading  the  system  would  again  be 
activated  to  maintain  the  desired  heading. 

A  block  diagram  of  the  data  collection  and  orient- 
ation systems  is  shown  in  Figure  2. 

Preliminary  data  were  obtained  during  February, 
March,  and  April  of  1966  while  Flip  was  operating  in 
Hawaiian  waters.   Analysis  of  these  data  indicated  the  need 
to: 

1.  increase  the  sensitivity  and  reliability  of  the  iso- 
therm followers 

2.  incorporate  the  thermometer  staffs  as  verniers 

3.  change  the  data  recording  system  to  incorporate  2 

4.  reposition  the  booms  and  winches 

5.  change  the  cable  fairleads  because  of  4 

6.  replace  lightweight  telescoping  booms  with  heavy-duty 
ones 

7.  employ  the  orientation  system  to  keep  the  array  stabil- 
ized in  azimuth  for  as  long  a  time  as  possible. 

The  implementation  of  these  modifications  resulted 
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In  the  system  described  above.   A  short  seatrial  was  con- 
ducted near  Gatalina  Island  in  November.   Then  FLIP  was 
towed  to  deep  water  off  Baja  California  and  data  were 
collected  during  December  1966. 
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CHAPTER  III 
OBSERVATIONS 

A  station  southwest  of  San  Diego  (near  30°N,121°W) 
was  maintained  from  1  December  1966  to  14  December  1966. 
FLIP  was  allowed  to  drift  with  the  mean  integrated  current 
over  her  300  foot  draft.   Since  the  waves  travel  relative 
to  the  current,  this  drift  minimized  any  Doppler  effects. 
A  track  of  FLIP's  drift  is  shown  in  Figure  3. 

Due  to  the  aerodynamic  characteristics  of  FLIP's 
upper  structure,  she  tends  to  orient  with  her  keel  into  the 
wind.   Once  she  has  self-oriented,  the  propeller  system  has 
little  trouble  maintaining  the  heading.   If  the  wind 
direction  shifts  more  than  twenty  degrees,  however,  the  load 
becomes  too  great  on  the  propeller  system  and  re-orientation 
is  required.   It  was  possible,  nevertheless,  to  break  the 
data  into  15  short  series  during  which  the  variance  of 
array  orientation  was  small.   The  duration,  mean  azimuth  of 
the  keel,  and  azimuthal  standard  deviation  are  given  in 
Table  I. 

Time  records  of  series  4,  5,  and  6  are  presented  in 
Figure  4,  and  are  typical  examples  of  the  internal  oscil- 
lations during  the  entire  sampling  period.   Each  series 
contains  isotherm  depth  vs.  time  from  all  three  followers. 
The  middle  trace  corresponds  to  the  depth  scale  on  the  left. 
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Fisure  3, 


Chart  of  FLIP" 3  track  and  location 
of  generating  areas  (after  H.  W.  Menard, 
I960).   Numbers  on  track  refer  to  the 
day  in  December  1966  (0800  position). 
Letters  refer  to  generating  areas. 
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TABLE  I 

Division  of  Data  into  15  Series 

STANDARD 
SERIES  NO,   TIME  (Dec.  1966)      MEAN  HEADING  (°T)  DEVIATION ( ° T ) 


1 

061400-021120 

329.2 

4.8 

2 

021120-031005 

239.8 

3.2 

3 

021827-031712 

240.2 

3.1 

4 

031816-041701 

240.7 

2.6 

5 

041701-051547 

240.6 

1.2 

6 

051819-061704 

239.1 

2.8 

7 

070115-071633 

348.5 

2.7 

8 

071633-080625 

328.5 

4.1 

9 

080751-090636 

343.2 

4.3 

10 

090636-100515 

342.4 

1.9 

11 

100557-110443 

011.8 

2.3 

12 

101430-111315 

010.4 

3.5 

13 

120338-122108 

002.4 

27.3 

14 

122108-132003 

356.5 

6.8 

15 

132003-140556 

007.6 

10.6 

15 


ISU313W]    WU3H10SI    JO  Hld3C 


(SU3I3U1    HU3H10SI    JO  Hid30 


ISU3I3UI    UU3H:DSI    JO  Hld30 


Figure  4.   Depth  of  the  observed  iaotherm  vs  time. 

The  upper  and  lower  traces  have  been  offset  +  10  meters. 
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The  upper  and  lower  traces  have  been  offset  +  10  meters 
respectively.   The  interpolation  scheme  used  to  calculate 
the  isotherm  depth  may  be  found  in  Appendix  B. 

The  first  characteristic  one  should  notice  in  the 
time  records  is  the  predominance  of  the  tidal  oscillations. 
Tidal  frequencies  contain  most  of  the  energy  and  are,  there- 
fore, responsible  for  the  largest  vertical  motions  of  the 
isotherm. 

Superimposed  on  the  tidal  motions  are  the  high  fre- 
quency fluctuations  with  which  this  thesis  is  concerned. 
Several  interesting  properties  of  the  high  frequency 
motions  are  apparent  in  the  time  records.   First,  the 
records  from  the  three  elements  exhibit  good  visual  coher- 
ence throughout  the  data.   Second,  the  larger  amplitudes 
of  the  high  frequency  oscillations  seem  to  come  and  go  as 
if  they  were  modulated.   And  finally,  the  oscillations 
appear  to  be  always  present,  indicating,  at  least  in  a 
qualitative  manner,  that  they  may  be  stationary.   These 
conclusions  are  Justified  quantitatively  in  Chapter  V. 

The  second  property,  i.e.  that  of  amplitude  modu- 
lation, has  been  observed  in  shallow  water  by  LaFond  (1961) 
and  may  be  caused  by  two  effects.   The  most  likely  cause  of 
this  amplitude  fluctuation  is  the  beating  effect.   When 
energy  is  present  as  a  continuous  spectrum,  beating  occurs 
among  the  various  frequencies  and  a  rapidly  varying  ampli- 
tude envelope  is  developed. 
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The  second  possible  cause  is  a  resonance  condition 
described  by  Eckart  (1961)  and  hereafter  referred  to  as 
Eckart  Resonance.   When  the  waves  are  propagated  in  a  well- 
localized  group  —  localized  in  space,  frequency,  and  wave- 
number  —  a  resonant  coupling  between  the  layers  of  high 
stability  can  occur.   The  effect  of  this  resonance  is  to 
cause  the  energy  to  oscillate  between  the  layers.   If  the 
conditions  are  variable  along  the  propagation  path,  and  one 
measures  the  motion  in  one  layer  only,  an  amplitude  modu- 
lation will  be  observed,  with  the  modulating  frequency 
dependent  upon  the  sharpness  of  the  layers  and  their  verti- 
cal separation.   The  same  effect  would  be  observed  if  the 
source  is  not  stationary  in  a  spatial  sense.   That  the 
necessary  conditions  for  this  resonance  existed  is  shown  in 
Chapter  IV. 


CHAPTER  IV 
ANALYSIS  TECHNIQUES 

The  power  spectrum  has  become  a  basic  tool  for  pre- 
senting the  properties  of  long  complicated  time  series  in  a 
condensed  form.   Here,  the  method  of  first  Fourier  trans- 
forming and  then  summing  over  the  squares  has  been  employed. 

Consider  a  discretely-sampled,  real  series  of  finite 
length,  X(J),  J  =  0,  N  -  1.   The  X(j)  may  be  represented  by 
the  complex  series 

X(j)=Z*(m)ei27rim/N'    J-O.N-I  (i) 

m  =  o 

Then  the  complex  Fourier  coefficients  are: 

^-i!'x(i)ei2Tim%=o,N-i 

j  =  0 

Taking  into  account  the  symmetry  of  the  transform  of  a  real 

series,  a  spectral  estimate  of  elementary  frequency  bands 
with  two  degrees  of  freedom  is  obtained  by 

A     A*_x    _   ^    N 

2  (3) 


Sx(m)=  AmX(m)  X(m)  ,    m  =  0,^-l 


where  *  denotes  complex  conjunction  and  k^   =  1,  m=o 

Am.  =  2,  m>o 
To  increase  spectral  stability,  a  sum  over  adjacent  ele- 
mentary bands  is  formed. 
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Sx(hc)  =  Z  Sx(m)=  £Am$r(m)  *(m)  <4> 

m=p 

where  p  and  q  are  the  lower  and  upper  band  indices  over 

which  it  is  desired  to  sum  and  h        =  p  +  q  .   Letting 

center      p 

a  =  q  -  p  +  1  (the  number  of  bands  Included  in  the  sum), 
equation  ( 4  )  yields  2a  degrees  of  freedom  per  estimate  . 
Since  it  is  desirable  to  model  a  continuous  spectrum  of 

infinite  extent  which  has  dimension  length  /cycle/time  = 

2 

L  T,  it  is  necessary  to  divide  S(h)  by  the  band  width,  a. 

Thus 

Ex(hc)=%ij2}=¥^  AmX*(m)*(m)  (5) 

x   c       VJ        U  m=p 

and  simple  summing  is  replaced  by  band  averaging. 
The  conversion  from  the  indices  m,  hc,  a,  to  fre- 
quency is  accomplished  by  using  the  sampling  interval  A  t. 
Let  NAt  =  T 

then  azinu  2 7r(j At)  ir 

N  h  a 

so  that    f-GL  ,    fcenter=f  .and  Af  =  T 

A  further  discussion  of  these  techniques  is  pre- 
sented by  Tukey  in  Harris  (1966).   Haubrich  (1965)  shows 
that  these  methods  produce  the  same  result  as  those  of 
transforming  the  correlation  function.   The  actual  calcu- 
lations were  performed  on  a  GDC  3600  digital  computer  using 
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the  fast  Fourier  transform  techniques  of  Cooley  and  Tukey 
(1965). 

Once  the  series  has  been  transformed,  the  cross 
spectral  estimate  between  two  series  is  simply 

Gxy(hc)=^Z  AmX#(m)^(m)  <« 

7      u  m=p 

The  estimates  of  coherence  and  phase  difference 
are  respectively: 


Rxy(hc)  = 


Gxy(hc) 


Ex(hc)Ey(hc) 


(7) 


and 

</)xy  (hc)=  Arg(Gxy  (hc))y  y  leads  x 


Considering  the  triangular  array 

3 


X  ^(hc)=Z  Arg(G/hc)),;  =  l2,23,3 


(8) 


(9) 


For  a  single  train  of  plane  waves,  equation  (7)  is 
identically  equal  to  1  ,  and  equation  (9)  is  identically 
equal  to  0.   Loss  of  coherence  and  a  non-zero  2-»'/f  mav  De 
caused  by  a  finite  beam  width  (Munk  et  al,  1963;  Cox,  1962). 
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Cox  also  discusses  reduction  of  coherence  caused  by  turbu- 
lence and  variable  phase  velocity.   It  should  be  noted  that 
if  a  r  1  (i.e.  p  =  q)  equation  (7)  has  only  one  degree  of 
freedom  and  njlHc)  is  constrained  to  be  equal  to  one.   Sim- 
ilarly, 2^  '^(hc)  in  equation  (9)  is  constrained  to  be  equal 
to  zero.   Therefore  it  is  necessary  to  band  average. 
Failure  to  average  over  at  least  two  bands  would  lead  to 
the  conclusion  of  plane  wave  propagation,  when  some  other 
form  of  propagation  may  exist. 

Using  the  three  phase  differences  defined  in 
equation  (8),  it  is  possible  to  determine  three  vector  wave 

numbers.   Let  vector  wave  number  k  have  components  (k  ,  k  ). 

x   y 

The  three  phase  differences,  <jP^>   for  plane  wave  propagation 
become 

<Pj=   kX  >Q  +  ky  Xf  (10) 

where  X»  and  Yj  are  the  spatial  separation  components  be- 
tween elements  of  the  array.   Since  there  are  three 
equations  and  only  two  unknowns,  (kY,  kv),  the  wave  number 
is  over  determined.   Considering  the  equations  two  at  a 
time,  three  separate  solutions  for  the  vector  components 
may  be  found. 
Then  r  ,  ^       |2-|,/2 

krms=  1.3  L  M  J  (11) 


where 


km  =  ikx  +  Jky    m  =  1,3  (12) 
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i  and  J  being  unit  vectors  in  the  x  and  y  directions  and  m 
indexing  the  three  sets  of  solutions. 

The  mean  direction  of  propagation  is  defined  as 

0=  Arg  (k,  +  k2+~k3) 


(13) 


Beam  width  indication,^  9  ,  is  defined  as 


A0= 


max(Arg(ki))-min(Arg(kj)    i,j=l,3;i#j       a*) 


where  the  branch  cut  is  taken  to  form^e  as  the  minimum 
angle  between  vector  solutions  k,  and  k.  (as  in  figure  5). 

k, 


Figure  5 

Definition  of  A0 

It  should  be  noted  that  if  Rg  =  1  and,  thus, £"((>.  =  0, 
then  only  two  of  the  equations  are  independent,  the  k's 
are  identical,  and«^9  =  0.   These  results  agree  with  the 
earlier  discussion  of  plane  wave  propagation. 
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Once  a  vector  wave  number  with  a  small  A0  has  been 
found,  it  becomes  necessary  to  determine  whether  the  motion 
Is  due  to  internal  waves  or  to  some  other  mechanism  such  as 
anisotropic  turbulence.   The  distinction  between  wave  motion 
and  motions  due  to  other  causes  can  be  made  with  the  help  of 
diagnostic  dispersion  curves.   The  solutions, krm8   are, 
plotted  on  an  CU-K  diagram  and  tested  against  the  theoreti- 
cal dispersion  curves  to  see  if  there  is  agreement.   It  must 
be  remembered,  however,  that  solutions  which  fall  the  test 
may  still  be  wave  motions  not  due  to  single  plane  waves  (i.e. 
multiple  modes,  or  two  or  more  plane  waves  of  the  same  mode 
and  frequency  traveling  in  different  directions,  etc.). 

To  obtain  the  dispersion  curves,  it  is  first  neces- 
sary to  calculate  density,   p   (z),  from  the  observed 
temperature,  T(z),  and  salinity,  S(z),  distributions.   Next 
p\Z)   is  used  to  calculate  N\Z/  .   Finally,  for  each  fre- 
quency band  of  interest,  the  wave  equation  is  numerically 
integrated  to  solve  for  the  kn. 

Figures  6  and  7  show  the  observed  temperature  and 
salinity  distributions,  the  calculated  density  and  stability 
frequency  distributions,  and  the  theoretical  dispersion 
curves  for  the  first  10  modes.   Figure  6  reflects  conditions 
Just  prior  to  the  commencement  of  wave  data  collection,  and 
Figure  7  reflects  the  conditions  immediately  after  termi- 
nating data  collection.   The  temperature  distribution  in  the 
upper  200  m  was  obtained  by  manually  lowering  an  isotherm 
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Figure  6.   Propagation  conditions  for  1  December  1966. 
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Figure  7.   Propagation  conditions  for  14  December  1966. 
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follower.   The  distributions  of  salinity  and  of  temperature 
below  200  m  were  obtained  from  nearby  CalCOFI  stations 
(Figure  6  in  December,  Figure  7  in  May).   The  difference  in 
the  two  salinity  distributions  at  shallow  depths  is  probably 
attributable  to  sampling  at  different  seasons  as  opposed  to 
being  a  permanent  spatial  feature.   In  any  case,  temperature 
is  the  controlling  factor  within  the  thermocline  and  the 
salinity  distributions  were  used  for  "order  of  magnitude" 
purposes.   Since  spatial  aliasing  will  occur  for  X  >  .016, 
the  dispersion  curves  are  terminated  at  this  value.   Below 
400  meters  T(z)  and  N(z)  monotonlcally  decrease  while  S(z) 
and  p  (z)  monotonlcally  increase  with  increasing  depth.   The 
values  of  these  parameters  corresponding  to  figures  6  and  7 
are  tabulated  in  Table  II. 

Several  features  of  Figures  6  and  7  are  worth 
further  consideration.   Both  stability  frequency  distri- 
butions have  many  relative  maxima.   These  maxima  are 
reflected  in  the  dispersion  curves  as  wiggles  (Fig.  6,  3rd. 
mode)  or  changes  in  curvature  (Fig  6,  2nd  mode).   It  is  the 
coupling  between  these  layers  of  high  stability  which  pro- 
duces Eckart  Resonance.   The  greater  the  vertical  ■ 
separation  between  the  layers,  the  sharper  the  resonance. 
The  resonance  regions  are  shown  in  the  dispersion  curves  as 
areas  where  two  modal  lines  converge  and  then  diverge.   For 
example,  in  Figure  7,  a  resonance  region  between  the  second 
and  third  mode  is  apparent  at  a  frequency  slightly  greater 


than  4cph  (X-1  =  .0085  )•  The  coupling  seems  to  occur  be- 
tween the  layers  centered  about  75  meters  and  250  meters. 
Other  regions  of  resonance  are  also  present. 

TABLE  II 
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Distribution  of  Temperature,  Salinity,  Density, 

and  Stability  Frequency  below  400. 

1  December  1966 

Depth  Temperature   Salinity  Density  Stability  Frequency 

Ccph) 
2.88  ^^ 
2.28 
2.06 
1.93 
1.77 
1.76 


500 

5.93 

34.20 

1.02993 

1000 

3.91 

34.47 

1.03320 

1500 

2.76 

34.57 

1.03622 

2000 

2.10 

34.63 

1.03912 

3000 

1.63 

34.66 

1.04472 

4000 

1.56 

34.66 

1.05021 

14  December  1966 


500 

5.92 

34.22 

1.02990 

3.00 

1000 

3.96 

34.49 

1.03321 

2.28 

1500 

2.89 

34.58 

1.03622 

2.06 

2000 

2.12 

34.64 

1.03913 

1.93 

3000 

1.65 

34.68 

1.04474 

1.80 

4000 

1.58 

34.69 

1.05024 

1.77 

The  two  sets  of  stability  frequency  and  dispersion 
curves  show  several  interesting  differences.   The  upper 
high  stability  layer  is  sharper  and  slightly  deeper  in 
Figure  7  relative  to  figure  6.   On  the  other  hand,  the 
lower  layer  is  weaker  and  much  deeper,  increasing  the 
separation  between  layers  and  causing  the  2nd  -  3rd  mode 
resonant  region  to  shift  to  a  lower  frequency  and  inverse 
wavenumber.   While  the  first  modal  curve  shifted  only 
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slightly  and  only  in  the  low  frequency  region,  all  of  the 
higher  modal  curves  shifted  toward  lower  frequencies.   The 
weakening  of  the  intermediate  relative  maxima  may  also  ex- 
plain why  the  curves  of  Figure  7  are  much  smoother  than 
those  of  Figure  6.   Notice,  however,  that  the  regions  of 
slightly  negative  curvature  still  remain  in  modes  2  and  3. 

These  differences  may  have  been  caused  by  both 
temporal  and  spatial  changes  in  the  internal  structure. 
But,  because  the  internal  waves  are  continually  changing 
the  subsurface  structure,  it  is  possible  that  the  mean  con- 
ditions were  stationary  throughout  the  data  collection 
interval »  and  that  the  differences  between  figures  6  and 
7  were  caused  by  taking  the  two  profiles  when  different 
phase  relationships  existed  among  the  waves.   In  this  case 
these  figures  would  represent  perturbations  about  the  mean 
conditions  at  the  two  spatially  and  temporally  separated 
points. 

Figure  8  shows  the  mean  depth  of  the  16.85°C  iso- 
therm for  series  1-12  and  the  16.00°C  Isotherm  for  series 
13-15  at  each  of  the  three  sea-sensing  elements.   The 
16.85°C  Isotherm  shows  large  oscillations  during  the  first 
12  series.   For  the  last  three  series,  however,  it  was 
necessary  to  shift  to  a  lower  temperature  in  order  to  re?- 
main  within  the  thermocline.   Therefore,  it  is  probable  that 
the  differences  are  due  to  a  combination  of  both  mechanisms. 
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CHAPTER  V 
RESULTS  OF  ANALYSIS 

The  fifteen  series  tabulated  In  Table  I  were 
analyzed  by  the  methods  of  Chapter  IV.   As  examples,  the 
results  for  series  1,  4,  8,  10,  and  14  are  shown  graphically 
In  Figures  10  through  14.   The  layout  of  each  figure  la  as 
follows; 

The  upper  left  hand  graph  shows  the  three  energy 
spectra  -  one  for  each  Isotherm  follower.   Spectral  averag- 
ing (equation  (5)  )  was  done  over  eight  elementary 
frequency  bands  yielding  16  degrees  of  freedom  per  spectral 

estimate.   For  16  degrees  of  freedom,  there  is  80$  confi- 

2      2  2  2      2 

denoe  that    ,58o  <S  <1.47cr     and  95%  confidence  that   .43cr  *S 

2        2  ? 

<1.80cf  where  S  is  the  spectral  estimate  and  d     is  the 

true  spectral  value.   Continuing  down  the  left  hand  side  of 
the  page,  the  second  graph  shows  the  three  pair  coherences 
as  defined  in  equation  (7).   For  16  degrees  of  freedom,  the 
95$  significance  level  for  coherence  is  .25.   The  third 
graphs  display  5  $   as  defined  in  equation  (9).   The  upper 
right  hand  graph  shows  the  computed  wave  number  magnitudes 
(equation  (11)  )  plotted  as  points  and  the  theoretical  dis- 
persion curves  (modes  1  to  10)  of  1  December  1966  plotted 
as  lines.   Since  the  modal  lines  converge  in  the  lower  fre- 
quency region,  an  expanded  scale  plot  for  series  1  is  shown 
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Figure   9.      Expanded   scale  dispersion  diagram   for  series   1. 
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FIGURE  10     SERIES  I 

ENERGY  SPECTRA*  PAIR  COHERENCES* 
AND 
SINGLE  WAVE  SOLUTIONS 
011400-021120  DECEMBER  1966 
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FIGURE  II      SERIES  4 

ENERGY  SPECTRA?  PAIR  COHERENCES? 
AND 
SINGLE  NAVE  SOLUTIONS 
031816-041701  DECEMBER  1966 
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FIGURE  12     SERIES   8 

ENERGY   SPECTRA?    PAIR  COHERENCES? 
AND 
SINGLE  WAVE  SOLUTIONS 
071633-080625  DECEMBER   1966 
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FIGURE  13      SERIES   10 

ENERGY  SPECTRA?    PAIR  COHERENCES? 
AND 
SINGLE  WAVE  SOLUTIONS 
090636-100515  DECEMBER   1966 
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FIGURE  14      SERIES   14 

ENERGY  SPECTRA?  PAIR  COHERENCES? 
AND 
SINGLE  WAVE  SOLUTIONS 
122108-131954  DECEMBER  1966 
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in  Figure  9  (corresponds  to  boxed  area  of  figure  10)  .   The 
final  graph  displays  the  directional  information.   Each  line 
runs  from  (6--^|)  to  (e+ — §)  where  a  and  A6  are  defined  by 
equations  (13)  and  (14)  respectively. 

Various  narrow  band  spectral  peaks  are  apparent 
throughout  the  15  series.   However,  the  peak  heights  lie 
between  the  80$  and  95$  confidence  limits  and  the  fre- 
quencies at  which  they  occur  change  from  series  to  series. 
The  peaks  are  caused  by  nearby  frequencies  having  the 
proper  phase  relationships  to  accentuate  their  particular 
frequency  band  during  the  short  duration  of  sampling.   This 
result  is  consistent  with  the  expected  results  from  sampling 
a  continuous,  random,  Gaussian  process. 

The  coherence  is  quite  high  for  the  lower  frequen- 
cies but  in  some  cases  falls  off  quite  rapidly  in  the  5  to 
7  cph  region.   Similarly, £  $   remains  close  to  zero  when  the 
coherence  is  high  and  then  appears  to  oscillate  randomly 
as  the  coherence  drops.   The  rate  at  which  the  coherence 
falls  is  greatest  when  the  energy  spectrum  is  also  falling 
rapidly  with  frequency  and  vice  versa.   A  possible  expla- 
nation for  the  fall  off  is  that  in  these  cases  there 
appears  to  be  a  lack  of  first  mode  solutions  at  the  higher 
frequencies.   For  modes  greater  than  1,  the  wavelength 
becomes  smaller  than  the  array  spacing  near  these  frequen- 
cies.  Therefore  as  frequency  increases  more  and  more 
wavelengths  occur  between  the  elements  of  the  array  (see 
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Cox  1962).   In  other  cases,  however,  the  coherence  rises 
again  to  quite  high  values  and  then  falls  off  in  the  10  cph 
region  with  £  <f>  behaving  in  a  consistent  manner.   In  these 
cases,  there  appears  to  be  a  predominance  of  first  mode 
solutions  above  5  cph  (compare  series  1  and  8). 

Some  of  these  solutions,  particularly  above  5  cph, 
lie  to  the  right  of  the  first  mode  dispersion  curve  but 
have  good  coherence,  small  J.  <$,  and  narrow  beam  width,  AG. 
In  several  series,  these  solutions  actually  seem  to  form  a 
separate  modal  line  below  that  of  the  first  mode.   The  best 
example  can  be  found  in  series  8  where  there  is  a  multitude 
of  narrow  beamed  high  frequency  waves  traveling  toward  270°. 
The  most  likely  cause  of  this  phenomenon  is  a  Doppler  shift 
due  to  current  shear.   FLIP  drifts  with  the  mean  Integrated 
current  over  her  entire  300  foot  draft  and  these  high  fre- 
quency waves  are  superimposed  on  the  mean  integrated  current 
over  their  relatively  narrow  high  stability  layer  (see  Fig- 
ures 6  and  7).   The  difference  between  these  currents 
required  to  account  for  the  observed  Doppler  shifts  is  on 
the  order  of  only  5  cm/sec.   The  details  of  the  necessary 
calculations  and  a  specific  example  may  be  found  in 
Appendix  C. 

The  solutions  considered  acceptable  as  single  wave 
propagation,  together  with  their  modal  numbers  and  mean 
propagation  directions  are  listed  in  Table  III.   Solutions 
falling  below  mode  1  are  included  as  mode  1  and  those 
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TABLE    III 

Summary   of  Acceptable   Solutions 

Band  No.  2        3  4  5  6  7  8  9        10       11  12        13        14        15       16       17        18        19       20       21        22       23        24       25       26        27        28        29         30  31  32  33  34         35  36  37 

FREQUENCE 

CPH .33    .68   1.03   1.38   1.74   2.09  2.44   2.79  3.14  3.49     3.94  4.20  4.55  4.90   5.25  5.60   5.96  6.31   6.66  7.01   7.36  7.71  8.06  8.42  8.77   9.12    9.47  9.82  10.17  10.52   10.88  11.23  11.5?   11.93     12.28     12.63 
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falling  above  mode  10  are  Included  as  mode  10.   All  other 
modal  Indexes  refer  to  the  nearest  modal  line.   The 
solutions  were  also  tested  against  the  14  December  dis- 
persion curves.   Complete  modal  agreement  was  found  except 
at  4.9  cph  in  series  8.   Here,  the  1  December  curves  indi- 
cate third  mode  and  14  December  curves  indicate  second  mode 
but  at  a  2nd-3rd  mode  resonant  region. 

Although  high  modes  frequently  occur  at  the  low 
frequencies,  the  first  mode  predominates  in  the  high  fre- 
quency range.   This  feature  is  probably  caused  by  the  wave 
guide  narrowing  with  Increasing  frequency  -  thus  subjecting 
the  higher  modes  to  greater  shear  until  they  eventually 
become  unstable.   In  general,  it  is  felt  that  the  modal 
solutions  are  quite  good.   Notice  the  good  agreement  between 
the  wave  number  solutions  and  the  theoretical  modal  lines  in 
Figure  9.   Credence  to  this  belief  is  also  given  by  the  good 
fit  obtained  in  the  slowness  plots  of  the  next  chapter.   It 
must  be  remembered  however  that  any  particular  modal  so- 
lution, and  especially  those  with  large  wave  numbers,  may 
have  been  shifted  relative  to  the  dispersion  curves  by  the 
Doppler  effect  discussed  above.   Table  IV  lists  the  number 
of  solutions  in  each  mode. 
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TABLE  IV 

Number  of  Acceptable  Solutions  In  Bach  Mode 

Percent  of  Solutions 
66 
12 
6 
3 
2 
1 

1/2 
1/2 
1 
8 


100# 


Numerous  wave  solutions  which  appear  to  be  station- 
ary over  several  days  may  be  found  by  referring  to  Table  III 
For  example  the  solutions  for  fc  =  2.44  cph  show  a  2nd  mode 
Internal  wave  traveling  153  -  1°  during  series  1,  2,  and  3. 
Similarly,  for  f  =  7.71,  the  solutions  for  series  5,  6,  7, 
8,  and  9  show  9  =  265  -  15°  and  for  tQ  =  8.06,  e=  279  •  3° 
during  series  8,  9,  and  10. 

To  define  the  general  shape  of  the  energy  spectrum, 
the  45  spectra  (three  for  each  series)  were  ensemble 


Mode 

Number  of  Solutions 

1 

257 

2 

47 

3 

23 
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11 
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9 

6 
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7 

2 

8 

2 

9 

5 

10 

_2i_ 

Total 

392 
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averaged.   Considering  that  the  three  spectra  from 
each  series  are  not  independent,  and  allowing  for  the  over- 
lapping in  time  of  several  of  the  series,  the  resulting 
spectral  estimates  have  194  degrees  of  freedom.   For  194 

degrees  of  freedom  there  is  80$  confidence  that  ,87o2  <S2 

2  2  2      2 

*1.13o*  and  95%  confidence  that  .8lcr  *S  <1.21o*  .   A  log-log 

plot  of  the  average  spectrum  is  shown  in  Figure  15. 

The  spectrum  has  been  broken  Into  characteristic 

regions  according  to  slope.   The  slope  of  the  curve  and  the 

power  of  frequency  at  which  the  spectrum  falls  are  Identical. 

Initially  (.2  cph  to  1.7  cph)  where  the  wave  guide  (W<|\)(z)) 

occupies  the  entire  water  column,  the  spectrum  falls  as 

-2  2  i 

f   *  .   At  1.7  cph,  the  slope  changes  to  f  "*-L  and  the  wave 

guide  begins  to  narrow  with  succeeding  spectral  estimates 
(i.e.,  for  *c   =  1.74  cph,  Az  =  4000m   for  fc  =  2.09  cph, 
Az  =  1430  m,  and  for  fc  =  2.44  cph.  Az  =860m;  where  A  z  is 
the  thickness  of  the  wave  guide).   From  5  cph  to  12  cph, 
the  spectrum  falls  as  f  ~3«25#   Above  12  cph,  no  Internal 
waves  appear  to  be  present   (the  4  acceptable  solutions 
above  12  cph  in  Table  III  have  been  Doppler  shifted  ),  and 
the  rate  of  fall  of  the  spectrum  begins  to  decrease  until 
at  50  cph  it  reaches  the  Inherent  noise  level  of  measure- 
ment.  This  region  of  zero  slope  indicates  a  root  mean 
square  noise  amplitude,  £  rms  (hc),  of  .014  meters. 

The  height  of  the  peak  at  4.2  cph  lies  between  the 
80$  and  95$  confidence  limits.   Also, the  individual  series 


43 


10 


>-    10 

C/) 


UJ 
Q 

cr 

UJ 
UJ 


10 


10-2. 


10 


-3 


T 1 — I     I    ■    I  I 1 1 1 1 — |     |    |    |  | 1 1 1 1 — |     |    I   I 


-3.25 


I0"41    i    i    i  i  i  1 1 i       i     i i  i  i  i  i  i i i I    I   M  l  l  J 


I  10 

FREQUENCY  (CPH) 


100 


Figure  15.   Average  energy  spectrum. 
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solutions  for  this  frequency  do  not  indicate  any  apparent 
reason  for  the  peak's  existence.   Therefore,  it  is  felt 
that  the  peak  has  no  statistical  significance. 

The  region  between  1.7  and  5  cph  coincides  with  the 
region  of  good  wave  solutions  and  the  energy  appears  to  be 
attributable  to  Internal  waves.   Using  a  Richardson's 
number  instability  argument,  Phillips  (1966)  has  proposed  a 
theoretical  equilibrium  spectrum  for  first  mode  propagation 
within  the  thermocline  for  frequencies  such  that  the  stabil- 
ity layer  is  small  compared  to  the  total  depth.   His 
spectrum  is  proportional  to   f    for  kd «  1  and  to  f  "^ 
for  kd  »  1  (where  d  is  the  thermocline  depth).   Since  pre- 
ceding analysis  has  shown  many  solutions  of  mode  higher  than 
one  within  the  region  kd  «  1,  the  agreement  between  the  ob- 
served spectrum  and  the  equilibrium  spectrum  indicates  that 
the   f    power  law  may  apply  in  general  for  kd<<  1,  and  is 
not  limited  solely  to  the  first  mode.   The  change  in  slope 
from   f  "*  to   f-3.25  occurs  where  there  is  a  rapid  narrow- 
ing of  the  wave  guide,  and,  as  Phillips  predicted,  where 
kd  =  1  (see  Figures  6  and  7).   For  kd»l,  the  first  mode 
is  predominant  and  the  agreement  between  the  observed 
spectrum  (   f  "3.25)  and  the  equilibrium  spectrum  (  f  ~^) 
furnishes  additional  support  to  Phillips'  Theory. 


CHAPTER  VI 

DIRECTIONAL  PROPAGATION  AND  POSSIBLE  SOURCES 

To  determine  whether  there  were  any  preferred 
directions  of  propagation,  the  accepted  solutions  (Table 
III)  were  subjected  to  further  analysis.   Each  solution  was 
spread  out  over  its  beam  width  and  normalized  so  that 

-360° 

J    SWdxjJ   =  I  (15) 

where  SWO=I/A0  ;  0-A  9/2  <  +  <  9  +A  9/2 
S(i//)=  0   ;  ^<9-A9/2,y>9+A9/2 

The  solutions  were  then  quantized  in  one  degree  intervals 
and  added  together. 

P(J)-ZSj('j)  ;  j-l/2,  3/2,... 719/2        (16) 
i 

where  i  Indexes  the  solutions  and  J  represents  the  mid- 
direction  of  the  interval.   Finally  a  10  degree  wide  run- 
ning mean  filter  was  used  to  smooth  the  results. 

Q(f)-I  P(mU=l,2...360;  P(-n)-P(360-n)     (17) 

The  resulting  beam  patterns  are  shown  in  Figure  16.   The 
top  pattern  includes  all  fifteen  series.   The  remaining  5 
patterns  include  three  series  each  and  show  the  changes 
which  occur  with  time.   Several  strong  lobes  are  noticeable 
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Figure  16.      Horizontal   beam  patterns. 
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within  the  patterns.   The  prominent  lobe  toward  270°  pro- 
vides an  interesting  example.   It  is  virtually  absent  during 
the  first  3  series.   It  begins  to  appear  during  the  next  3 
series,  and  for  series  7  to  9,  overshadows  any  other 
direction.   It  then  decreases  but  remains  present  for  the 
final  6  series.   The  existence  of  these  lobes  leads  to  the 
conclusion  that  preferential  directions  indeed  do  exist. 
The  existence  of  directional  propagation  leads,  in 
turn,  to  the  consideration  of  possible  sources.   Azimuthal 
information  is  provided  by  the  beam  patterns.   To  determine 
range,  the  dispersive  character  of  the  wave  propagation  is 
utilized.   The  wave  energy  travels  at  the  group  velocity, 

n*  '   3k   ^Xn  (where  n  indexes  the  mode) 

which  is  the  slope  of  the  local  tangent  to  the  dispersion 
curves  at  (f,|/ An).if  the  source  is  localized,  so  that  the 
distance  of  travel  is  the  same  for  all  wave  groups,  then  for 
two  wave  groups  arriving  at  t  and  t  +4t  t: 

X  =  G(f0)t  (is) 

X=G(f,)(t+At)  d9) 

so  that       . 

XG"(fo)=t  (20) 

XG"'(f,)  =  t+At  (21) 

Letting  AG~  =G  (f|)-G  (fo) and  subtracting  equation  (20) 
from  equation  (21)     X  =  At/AG"'  (22) 
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Thus,  range  to  a  source  is  obtained  by  making  slowness  plots 
(time  of  arrival  vs.  inverse  group  velocity).   The  slope  and 
t- intercept  of  a  reasonable  straight  line  fit  give  range  and 
time  of  generation  respectively.   The  slowness  plots  for  the 
lobes  centered  at  270°,  070°  and  025°  are  shown  in  Figures 
17,  18  and  19.   The  slope  for  270°  (=  1.59  day-meters/sec 
gives  a  distance  of  approximately  75  miles.   Looking  in  the 
reciprocal  direction  (i.e.,  090°)  at  a  distance  of  75  miles 
from  FLIP's  track,  one  observes  a  prominent  ridge  of  sea- 
mounts  (labeled  "A"  in  Figure  3). 

In  Figure  18  several  different  slopes  are  apparent. 
The  steepest  slope  (7.14  day-meters/sec)  indicates  a  distance 
of  330  miles  and  again,  the  source  seems  to  be  a  cluster  of 
seamounts  (labeled  "BM  in  Figure  3).  The  next  steepest 
slope  (2.18  day-meters/sec)  indicates  a  distance  of  100 
miles  and  the  seamount  marked  "C"  in  Figure  3  seems  to  be 
the  source.   The  least  slope  (.92  day-meter/sec)  indicates 
a  distance  of  43  miles.   Since  this  solution  occurs  during 
the  time  when  FLIP  was  rapidly  drifting  toward  the  south,  it 
is  felt  that  the  seamount  at  MD"  (Figure  3)  was  the  source 
and  that  the  wave  group  also  was  subjected  to  advectlon  in  a 
southerly  direction. 

In  Figure  19,  only  the  slope  of  the  heavy  line 
(5.00  day-meters/sec)  corresponds  to  a  topographic  feature 
(range  =  233  miles  -  area  labeled  "E"  in  Figure  4).   The 
other  lines  give  ranges  of  110,  345,  and  24  miles.   They  are 
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Figure  18.   Slowness  plot  for  direction  070°.   Numbers  refer 
to  frequency  band  numbers  of  Table  III.   Super- 
scripts indicate  mode. 
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Figure  19.   Slowness  plot  for  direction  025  .   Numbers  refer 
to  frequency  band  numbers  of  Table  III.   Super- 
scripts Indicate  mode. 
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shown  to  indicate  that  many  of  the  straight  line  fits  give 
ranges  to  which  no  correlation  with  topography  was  found. 

The  quantity  of  solutions  used  to  determine  the 
slopes  makes  the  range  determinations  subject  to  doubt. 
These  results  were  presented  not  as  absolute  proof  that 
seamounts  act  as  internal  wave  sources,  but  merely  to 
demonstrate  that  the  solutions  point  to  them  as  likely  pos- 
sibilities.  This  conclusion  is  supported  by  Sandstrom 
(1966)  who  has  shown  that  topographic  features  tend  to 
increase  the  energy  in  internal  waves. 


CHAPTER  VII 
CONCLUSIONS 

High  frequency  internal  waves  in  the  deep  ocean  are 
characterized  as  a  broad  band  phenomenon  with  a  continuous 
distribution  in  frequency.   The  spectral  shape  is  consistent 
with  Phillips'  theory  of  a  shear  limiting  process,  and  the 
predominance  of  the  first  mode  at  high  frequency  gives 
additional  emphasis  to  the  importance  of  shear  instability 
in  Internal  wave  propagation. 

Over  short  intervals  of  time,  a  narrow  frequency 
band  within  the  continuum  displays  great  simplicity  in 
modal  and  directional  structure  and  can  be  adequately  de- 
scribed as  a  horizontally  plane  wave  (i.e.,  small  beam 
width)  of  a  single  vertical  mode.   In  addition,  it  has  been 
shown  that  this  simplified  structure  is  statistically 
stationary  for  time-spans  on  the  order  of  several  days. 
This  structural  simplicity  indioates  that  the  waves  are  well 
dispersed  and  implies  the  existence  of  fairly  localized 
(both  in  time  and  space)  sources. 

Considering  the  properties  of  the  total  (broad  band) 
wave  field,  additional  evidence  of  localized  sources  is 
found.   The  horizontal  nature  of  the  wave  field  displays 
definite  anlsotropy,  Indicating  well-defined  directions  of 
narrow-beamed  wave  propagation.   While  the  exact  character 
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of  the  anisotropy  may  change  appreciably  In  times  of  the 
order  of  several  days,  well  defined  directions  of  narrow- 
beamed  propagation  always  appear  to  be  present. 

The  time-history  of  the  wave  field  displays  the 
presence  of  dispersive  wave  trains.   As  well  as  indicating 
localized  sources,  their  presence  makes  possible  the  calcu- 
lation of  source  range.   These  range  determinations,  coupled 
with  the  well-defined  propagation  directions,  indicate  that 
nearby  topographic  features  such  as  seamounts  act  as  gener- 
ators of  high  frequency  internal  waves. 

In  view  of  the  wide  distribution  of  seamounts 
throughout  the  oceanic  basin,  the  superposition  of  arrivals 
from  the  multitude  of  sources  should  produce  a  confused 
"internal  sea".   Therefore,  the  waves  must  not  persist  over 
long  propagation  paths. 
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APPENDIX  A 


Theory  of  Small  Amplitude  Internal  Waves 
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I.   List  of  Symbols 

C  velocity  of  sound 

f  inert lal  frequency 

9  acceleration  of  gravity 

k  wave  number  magnitude 

Kx  x-component  of  wave  number  vector 

ky  y- component  of  wave  number  vector 

n  mode  index 

time 

X  horizontal  coordinate 

7  horizontal  coordinate  orthogonal  to  x 

Z  vertical  coordinate  (positive  up) 

D  depth  of  bottom 

F  force 

N  stability  frequency 

q  specific  volume 

£  vertical  displacement 

p  density 

OJ  angular  frequency 

Q V  Infinitesimal  volume 

AZ  small  vertical  displacement 

(  )  second  derivative  with  respect  to  time 

\   )-7j  second  derivative  with  respect  to  Z 
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II.   Infinites imal  Waves 

The  following  presentation  consists  of  only  those 
aspects  of  Internal  wave  theory  necessary  to  the  understand- 
ing of  the  observational  and  analytical  techniques  employed 
in  this  dissertation.   More  detailed  formulations  may  be 
found  in  Eckart  (I960)  and  Tolstoy  (1963). 

Consider  a  vertically  stratified,  compressible  fluid 
at  rest  in  a  gravity  field.   If  a  parcel  of  fluid  is  raised 
vertically  from  its  equilibrium  position  and  then  released, 
it  will  tend  to  sink  with  the  force: 

F=g(/H/0+!£Az)dV)7,^  AzdV  (ia) 

the  first  term  being  due  to  buoyancy  and  the  second  compres- 
sibility.  The  vertical  equation  of  motion  becomes: 

yodVz—g^zdV-^zdV  (2A) 

For  infinitesimal  motions,  the  solution  to  equation  (2)  is 
that  of  a  harmonic  oscillator  with  the  frequency: 


N<W-i£?-4 

Equation  (3A)  defines  the  Brunt-Vaisala  op  stability  fre- 
quency.  In  the  upper  ocean,  the  compressibility  term  is 
generally  small  compared  to  the  buoyancy  term  and  is  usually 
neglected.   Thus  we  have  harmonic  motion  caused  by  buoyancy 
forces  due  to  density  stratification. 
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Starting  with  the  three-dimensional  equations  of 
motion,  mass  conservation,  and  invoking  the  following 
assumptions: 

1)  N  =  N(z)  and  is  independent  of  horizontal 
coordinates 

2)  Mean  velocities  are  equal  to  zero 

3)  Boussenesq  Approximation  -  Density  variation 
may  be  neglected  except  where  it  causes  buoyant  forces 

4)  Nonlinear  terms  are  negligible 

we  arrive  at  the  first  order  equation  of  motion 

N2(z)-co2 


t>zz-U^-  v„C=  0 


(4A) 


where  Cl)  is  the  angular  frequency  and   is  the  inertial  fre- 
quency due  to  rotation  of  the  earth  and  £  is  the  vertical 
displacement  from  equilibrium. 

If  we  assume  a  solution  of  plane  horizontally  pro- 
gressive waves  of  the  form 

£(x,y,z,t)  =  £(z)e  y 


(5A) 


and  k*  =  k2  +  k2, 


then  (4A)  becomes 

2 


N  -or  1,2  1-  _  n 

C    +    2   X2  k  b  -  ^  (6A) 

Internal  waves  have  their  maximum  amplitude  below 
the  surface  of  the  water  column  and  at  this  maximum 
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£77<0 


zz 

but       ~*  £-L    =-K"    —^2      Tz 


2-u,2 


Therefore  either  f  2<0>2<N20r  N2<u;2<f  2.      m  the  ocean  N2»f2 
so   that  the  former  condition  is   the  correct  one,    and  f   and 
N  max  are  the  lowest  and  highest  frequencies  permitted  to 
free   internal  waves  of  form    (5A). 

If  CO   25>T    then  equation    (6A)   becomes 

Czz+K    cu2   t  -  0 

The  boundary  conditions  for  our  purposes  are 
£  =  0      at       Z  =  0       surface 

£  =  0     at      Z=~~D     bottom 
The  free  surface  boundary  condition  may  be  used  at  the  sur- 
face but  has  very  little  effect  on  the  internal  modes. 
For  arbitrary  N(z)  and  fixed  CO   ,  solutions  to 
equation  (7A)  are  found  by  iterative  numerical  techniques  - 
starting  with  C,   =  0  at  z  =-D  and  working  to  satisfy  the  sur- 
face boundary  condition.   For  a  given  CO     ,  solutions  can  be 
found  for  an  infinite  number  of  discrete  values  of  K.    For 

each  such  value  of  K  ,  there  is  a  corresponding  amplitude 

n 

function  C,   n.   The  index  n  relates  to  the  number  of  ampli- 
tude maxima  within  the  water  column,  and  (n-1)  to  the  number 
of  nodes  or  zero  crossings  between  the  maxima. 

Solutions  will  be  progressive  if  CO<  IMvZj somewhere  in 


the  water  column./"  (2)    will  be  "trigonometric"  where 
GJ<N(z)and  "exponential"  where  CJ>N(z)     (see  Figure  1A). 
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exponential 
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Figure  1A.   2nd  Mode  Amplitude  Function 

Within  the  thermocline  where  the  density  gradient, 
and  therefore  N(z),  is  a  maximum,  the  appropriate  vari- 
ation of  physical  parameters  and  their  effects  on  specific 
volume  are  given  in  Table  IA 


TABLE  IA 
The  Variations  of  Physical  Parameters  in  the 
thermocline  and  their  effects  on  specific  volume 

Parameter  Variation  X    a  £(  )           Specific  Volume  Change 

Temperature  -8°C  -2X10"  V°c               -1.6X10"3 

Salinity  +.3%>  -8X10"4/%o              -2.4X1CT4 
Pressure 


-4  A 

+100db  -4X10     /lOOdb        -4X10 
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It  has  already  been  shown  that  internal  waves  are  oscil- 
lations of  lsodenslty  surfaces.   Since  temperature  is  the 
dominant  parameter  within  the  thermocline,  the  motions  of  an 
lsodenslty  surface  may  be  approximated  by  those  of  an  iso- 
therm.  This  approximation  is  very  convenient  since 
temperature  is  much  easier  to  measure  than  density.   Temper- 
ature measurements  for  this  reason  have  been  widely  used  for 
obtaining  observational  data. 


APPENDIX  B 

Interpolation  Scheme 

The  following  parameters  were  sequentially  sampled 
from  each  sea-sensing  unit: 

Da»  Tl»  T2»  T3»  T4»  °b 
where  D  =  depth  and  T  =  temperature. 

Then  the  depth  D^  of  each  thermometer  when  sampled 

may  be  found  to  first  order  by 

ADi  =  (Db-Pa)i,  i  =  1,4 

5  (12) 

Di  =  Da  +  A^i  +  Si 

where  b.  is  the  vertical  offset  of  the  i   thermometer  from 
the  servo-system  thermistor. 

The  depth  of  the  isotherm  (DT)  to  first  order  is 
then  found  by  using  the  two  temperature  measurements  which 
bracket  the  isotherm: 

DT  =  Dt  +  [(D1+i  -  Bi)/(Ti.-T1+1)j  (TA  -  T) 

where  T  is  the  temperature  of  the  desired  isotherm. 
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APPENDIX  C 
Doppler  Effect 
Phase  velocity,  c,  is  defined  as 

c  =  w/k  =  f/X-' 

When  sampling  at  three  points  simultaneously,  the  Doppler 
shift  occurs  only  in  frequency,  leaving  the  wavelength 
unaffected.   Thus: 


qj  +  Aoj        f  +Af 
C  +AC  =  — r —  =  — ^T~ 


Ac  =   Af/X"1 

For  frequency  9.47  cph  of  series  8,  X  =8.62x10   ITI 
X-,=  8.62x  IO"3npToccurs  at  7.90  cph  on  the  first  modal 
curve  of  Figure  6. 


9.47-7.90   x    IP2      =     1.57    ^  5  Cm 
AC  =  :  8.62x  I0"3  "  3600  .31  Sec 
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